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METHOD AND SYSTEM FOR SIMULATING AN UNSTEADY FLOW FIELD 

TECHNICAL FIELD OF THE INVENTION 

This invention relates generally to the field of 
simulation systems and more specifically to a method and 
system for simulating a flow field. 
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BACKGROUND OF THE INVENTION 

Simulating the flow field around an aircraft 
transporting a store, such as a weapon, is important for 
the design of aircraft. Increasingly, modern aircraft 
5 transport their store internally within a closed bay. 

The bay is opened to deploy the store over a target. As 
the bay is opened, the flow field around the aircraft 
affects the deployment of the store, and may become 
unsteady around the open bay. Typically, the flow field 

10 is simulated to determine the effect of deployment of the 

store on the aircraft . 

Simulating the flow field around an aircraft, 
however, has posed challenges to aircraft designers. 
Wind tunnels may be used to simulate a flow field. Scale 

15 models of the aircraft are placed in a wind tunnel, and 

measurements of the flow field created in the wind tunnel 
are taken. The results from wind tunnel simulations, 
however, are often distorted since the size and mass of 
the scale model are different from that of the actual 

20 aircraft. Additionally, wind tunnel simulation is 

typically expensive. Computational simulations may also 
be used to determine the flow field. Calculating the 
flow field around an aircraft, especially an unsteady 
flow field, however, requires a large amount of computing 

25 power and time. Consequently, simulating the flow field 

around an aircraft has posed challenges to aircraft 
designers . 
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SUMMARY OF THE INVENTION 

While recent developments in aircraft simulation 
have provided improvements over prior approaches, the 
challenges in the field of simulation systems have 
5 continued to increase with demands for more and better 

techniques having greater effectiveness and efficiency. 
Therefore, a need has arisen for an improved method and 
system for simulating a flow field. 

In accordance with the present invention, a method 
10 and system for simulating a flow field are provided that 

substantially eliminate or reduce the disadvantages and 
problems associated with previously developed systems and 
methods . 

According to one embodiment of the present 
15 invention, a method for simulating flow field is 

disclosed. A grid comprising a plurality of cells is 
provided. Each cell is associated with a set of variables 
that describe a flow field. For each cell, a value for 
each variable is calculated from a previous value at each 
2 0 period for a predetermined number of periods. The 

calculated values for each variable are averaged to yield 
an averaged value for each variable. A flow field is 
determined from the averaged values. 

According to one embodiment of the present 

2 5 invention, a system for simulating a flow field is 

disclosed. A grid generator generates a grid comprising a 
plurality of cells. Each cell is associated with a set of 
variables for describing a flow field. A flow field 
module is coupled to the grid generator. The flow field 

3 0 module calculates a value for each variable from a 

previous value at each period for a predetermined number 
of periods. An averaging module is coupled to the flow 
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field module. The averaging module averages the 

calculated values for each variable to yield an averaged 
value for each variable. 

Embodiments of the invention may provide numerous 
5 technical advantages. A technical advantage of one 

embodiment is that values describing the flow field are 
averaged to smooth out instabilities, resulting in an 
accurate flow field simulation that uses relatively 
little computing power and time. Another technical 

10 advantage of one embodiment is that the values calculated 

using a time-warped process may be averaged in order to 
produce a simulation of sufficient accuracy, while using 
relatively little computing power and time. 

Other technical advantages are readily apparent to 

15 one skilled in the art from the following figures, 

descriptions, and claims. 
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BRIEF DESCRIPTION OF THE DRAWINGS 

For a more complete understanding of the present 
invention and for further features and advantages, 
reference is now made to the following description, taken 
5 in conjunction with the accompanying drawings, in which: 

FIGURE 1A illustrates an object for which a flow 
field may be simulated; 

FIGURE IB illustrates a grid used to simulate a flow 
field in accordance with the present invention; 
10 FIGURE 2 is a block diagram of a system for 

simulating a flow field in accordance with the present 
invention; and 

FIGURE 3 is a flowchart demonstrating a method for 
simulating a flow field in accordance with the present 
15 invention. 
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DETAILED DESCRIPTION OF THE DRAWINGS 

Embodiments of the present invention and its 
advantages are best understood by referring to FIGURES 1 
through 3 of the drawings, like numerals being used for 
5 like and corresponding parts of the various drawings. 

FIGURE 1A illustrates an example of an object for 
which a flow field may be simulated by the present 
invention. A flow field includes any fluid such as air 
or water, and an object includes any suitable solid in a 

10 flow field, for example, an aircraft 10 such as an F-22 

Raptor. Aircraft 10 includes a bay 12 used to transport 
and deploy a store 14 such as a weapon. Store 14, 
however, may include any object transported and deployed 
by aircraft 10, such as cargo. Bay 12 includes ejectors 

15 18 that eject store 14 over a desired target. Bay 12 

includes a door 16 that typically is closed during 
transportation of store 14 and opened to deploy store 14 . 
To simulate a flow field for aircraft 10, an electronic 
representation of aircraft 10 is generated using, for 

20 example, ACAD or other suitable solid modeling software 

program. 

FIGURE IB illustrates a grid 20 for aircraft 10 used 
to generate a flow field simulation. Grid 20 comprises 
cells 22 ci,...,c n that are illustrated as squares of grid 

25 20. Each cell 22 Ci is associated with a set of Q- 

variables Qi,...,Q m that describe the flow field at the cell 
22. "Each" as used in this document means each member of 
a set or each member of a subset of the set. Q-variables 
may include, for example, velocity variables such as an 

3 0 axial velocity variable, a transverse velocity variable, 

and a vertical velocity variable, a pressure variable, 
and a temperature variable. Q-variables, however, may 
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include any variables suitable for describing a flow 
field, such as momentum variables, a density variable, 
and an energy variable. Each variable Qk(t) takes on 
values q kt that describe the variable at a particular time 
5 t. For example, velocity values, a pressure value, and a 

temperature value describe the velocity, pressure, and 
temperature at a particular time. 

Grid 20 includes cells 22 of varying sizes. Larger 
cells may be used to describe regions with variables that 

10 change slowly, and smaller cells 22 may be used to 

describe regions with variables that change rapidly. For 
example, small cells 24 are used to describe the region 
inside bay 12 that has a complicated structure where the 
variables change relatively rapidly. Approximately 160 

15 cells 24 may span bay 12 in the axial direction, and 

approximately 45 cells 24 may span bay 12 in the lateral 
direction. Medium sized cells 26 are used to describe 
the regions just outside of bay 12, where variables 
change less rapidly. Large cells 28 are used to describe 

2 0 the regions far from bay 12, where variables change 

relatively slowly. 

The sizes of cells 22 are adjusted as the flow field 
changes. For example, grid 2 0 is generated for a closed 
bay 12, that is, bay 12 with door 16 closed, or for an 
25 open bay 12. Since a closed bay 12 has simpler geometry, 

grid 2 0 for closed bay 12 includes, for example, 
approximately 500,000 cells 22, whereas grid 20 for open 
bay 12 includes, for example, approximately 1,0 00,000 
cells 22. 

3 0 FIGURE 2 illustrates a system 10 0 for simulating a 

flow field for an object. System 100 may use any 
suitable configuration of components, software, and/or 
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hardware to simulate a flow field for an object. In the 
illustrated examples, system 100 includes an input-output 
(I/O) 110 that receives initial conditions for 
calculating a flow field, and outputs values for a flow 
5 field simulation. A processing unit 112 coupled to 

input-output 110 manages the process of and provides 
computational power for simulating a flow field. 
Processing unit 112 includes, for example, a combination 
of Hewlett-Packard V-class and X-class machines using 

10 eight to sixteen processors, or any other suitable 

combination of machines and processors. A memory 113 
coupled to processing unit 112 stores data during the 
simulation. A recorder 115 coupled to memory 113 records 
data in memory 113 in the form of reusable magnetic 

15 media. 

Software modules 108 perform calculations for 
simulating a flow field. A grid generator 114 coupled to 
processing unit 112 generates and adjusts grid 20 for 
aircraft 10. Grid 20 is generated with cells 22 c± 

20 suitably sized to sufficiently describe the flow field 

for aircraft 10. The sizes of cells 22 may also be 
adjusted as the flow field changes in time. Grid 
generator 114 includes computer aided design software 
such as PRO ENGINEER by PTC, INC. 

25 Flow field modules 115, such as a time-accurate 

module 116 and a time-warped module 118 coupled to grid 
adapter 114, calculate the values q kt for the variables 
Q k (t) in order to determine a flow field at time t. Time- 
accurate module 116 and time-warped module 118 calculate 

3 0 the values by applying flow field functions, such as 

Navier-Stokes flow field equations, to values in a time t 
to determine the values at a time t + 1. Flow field 
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modules 115 include computational fluid dynamics (CFD) 
software such as FLUENT, by FLUENT, INC. Values for the 
variables are determined at every time step, and the 
values are periodically recorded after a predetermined 
5 number of time steps. For example, values are calculated 

approximately every 2.0xl0" 5 seconds, and values are 
recorded at periods of every five time steps. The length 
of the time step is selected by balancing computation 
time with accuracy. A longer time step decreases 

10 computation time, but also decreases accuracy and 

numerical stability. If the desired time step can not be 
achieved during the initial calculation process, a 
shorter time step is used initially and then 
incrementally increased until it reaches the desired time 

15 step. Values are recorded as plot files stored in memory 

113. Plot files describe, for example, cell locations 
and values at the cell locations. 

Time-accurate module 116 determines the state of a 
flow field using a time-accurate process. Time-accurate 

2 0 module 116 calculates and records the values for all 

cells 22 at the same rate. Time-warped module 118 
determines the state of the flow field using a time- 
warped process. Time-warped module 118 determines and 
records the values for different cells 22 at different 
25 rates. For example, values of larger cells 28 are 

determined and recorded at a rate faster than the rate 
for values of smaller cells 24. 

Averaging module 12 0 coupled to time -accurate module 
116 and time-warped module 118 receives recorded values 

3 0 q kt from time-accurate module 116 or from time-warped 

module 118 and averages the values q kt for a predetermined 
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number p of periods to yield averaged values ~S ^ • 

For example, values q it for Q-variable Qi recorded for p 

1 p 

periods are averaged to yield averaged values ^ q H for 

Q-variable Qi . Averaging the values smooth out 

5 instabilities in the values, resulting in more accurate 

data . 

The number of periods is selected according to the 
accuracy of the recorded values and the desired accuracy. 
For example, averaging module 120 averages time-accurate 

10 values recorded for 200 to 400 periods. Averaging module 

12 0 may average time -warped values recorded for 
approximately 30 to 40 periods. Averaging module 120 
stores the averaged values in a plot file in memory 113. 
Averaging module 120 also gathers statistical information 

15 for the averaged values, such as standard deviation and 

root -mean- square values. 

Interpolation module 122 coupled to averaging module 
12 0 transforms the averaged values as associated with 
grid 2 0 to values associated with a survey grid, so that 

2 0 the averaged values may be compared with test data or may 

be used in other types of calculations. Software such as 
FIELD VIEW by INTELLIGENT LIGHT, INC., may be used to 
interpolate the values. Simulation software 124 coupled 
to interpolation module 122 uses the interpolated values 

2 5 in other simulation software. For example, store 

simulation software uses the values to determine how the 
flow field described by the values affects store 14 
transported by aircraft 10 upon ejection from aircraft 
10 . 
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FIGURE 3 is a flowchart illustrating a method for 
simulating a flow field. The method begins at step 212, 
where initial conditions for a flow field for aircraft 10 
are received. Initial conditions are calculated from 
5 user input describing, for example, the velocity of 

aircraft 10. The initial conditions may be selected to 
increase the speed of the calculation by estimating the 
final solution. For example, the values for the velocity 
variables are initialized to zero in the open bay 12. 

10 Grid 20 is generated for open bay 12 at step 214. Grid 

generator 114 adjusts cells 22 such that cells 22 are 
appropriate for describing a flow field around aircraft 
10. Grid 20 for open bay 12 may comprise, for example, 
approximately 1,000,000 cells 22. 

15 At step 216, the method determines whether a time- 

accurate process or a time-warped process is to be used. 
The time-accurate process is used if an accurate 
simulation is required. The time-accurate process, 

however, requires more processing time and power. The 

20 time-warped process produces a less accurate simulation, 

but also uses less processing power and time. For 
example, time-accurate module 116 may require the use of 
fourteen processors and approximately 14,000 CPU hours, 
while time-warped module 118 may require the use of only 

25 eight processors and approximately 115 CPU hours. 

If the time-accurate process is to be used at step 
216, the method proceeds to step 218 to initiate the 
time-accurate flow field calculation. The values q kt for 
the variables Qk are determined at each time step, and 

3 0 recorded periodically, for example, every five (5) time 

steps in a plot file. A transient period is determined 
at step 220. The transient period is determined by 
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examining the values. For example, the values of the 
pressure variable may spike during the transient period, 
and then become relatively steady after the transient 
period. The transient period may last for 1000 to 2000 
5 time steps, for example. 

After passing the transient period, the flow field 
calculation is continued at step 222. At step 224, the 
method determines whether a predetermined number of 
recording periods, for example, 200 to 400 periods have 

10 passed. If a predetermined number of periods have not 

passed, the method returns to step 222 to continue the 
flow field calculation. If the predetermined number of 
periods have passed, the method proceeds to step 226 
where the recorded values are averaged. For example, 

15 values from 2 00 to 4 00 periods are averaged. 

After averaging the recorded values, the method 
proceeds to step 22 7 to determine whether a predetermined 
final period has been reached. The final period signals 
the end of calculating the values. If the final period 

2 0 has not been reached, the method returns to step 222 to 

continue the flow field calculation. If the final period 
has been reached, the method proceeds to step 22 8 to 
interpolate the averaged values. Interpolation module 
122 adjusts the values as calculated for grid 2 0 to 

25 values associated with a survey grid in order to compare 

the results with test data at step 22 9. Other suitable 
operations may be performed on the interpolated values. 
For example, the values may be used in other simulation 
process . 

30 If a time-warped process is to be used at step 216, 

the method proceeds to step 230 to initiate the time- 
warped flow field calculation. The values q kt for the 
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variables Q k are determined at each step, and recorded at 
each period, for example, every fifty (50) time steps, in 
a plot file. A transient period is determined at step 
232 . The transient period is determined by animating the 
5 flow field from the plot files to illustrate the values. 

For example, the values spike during the transient period 
and then become steady after the transient period. The 
transient period may last for approximately 2000 time 
steps, for example. 

10 After passing the transient period, the flow field 

calculation is continued at step 234. At step 236, the 
method determines whether a predetermined number of 
recording periods, for example, 3 0 to 4 0 periods, have 
passed. If a predetermined number of periods have not 

15 passed, the method returns to step 234 to continue the 

flow field calculation. If the predetermined number of 
periods have passed, the method proceeds to step 23 8 
where the recorded values are averaged. For example, 
values from approximately 3 0 to 4 0 periods are averaged. 

2 0 After averaging the recorded values, the method 

proceeds to step 239 to determine whether a final period 
has been reached. If the final period has not been 
reached, the method returns to step 234 to continue the 
flow field calculation. If the final period has been 

2 5 reached, the method proceeds to step 24 0 to interpolate 

the averaged values. Interpolation module 122 adjusts 
the values as calculated for grid 2 0 to values associated 
with a survey grid in order to be able to use the values 
with other simulation software or to perform other 

3 0 operations on the results. 

At step 240, the interpolated values are used in 
other simulation software. For example, the values are 
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used in store simulation software in order to determine 
how the flow field described by the values affects store 
14 in bay 12 of aircraft 10. Other suitable operations 
may be performed on the values. For example, the values 
5 may be compared with test data. After using the values 

in the other simulation software, the method terminates. 

Although embodiments of the invention and its 
advantages are described in detail, a person skilled in 
the art could make various alterations, additions, and 
10 omissions without departing from the spirit and scope of 

the present invention as defined by the appended claims. 



